Double Binary Pulsar Test of Chern-Simons Modified Gravity 
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Chern-Simons modified gravity is a string-theory and loop-quantuni-gravity inspired effective 
theory that modifies General Relativity by adding a parity-violating Pontryagin density to the 
Einstein-Hilbert action multiplied by a coupling scalar. We strongly constrain non-dynamical Chern- 
Simons modified gravity with a timelike Chern-Simons scalar through observations of the double 
binary pulsar PSR J0737 — 3039A/B. We first calculate Chern-Simons corrections to the orbital 
evolution of binary systems. We find that the ratio of the correction to periastron precession to the 
general relativistic prediction scales quadratically with the semi-major axis and inversely with the 
square of the object's radius. Binary pulsar systems are thus ideal to test this theory, since periastron 
precession can be measured with sub-degree accuracies and the semi-major axis is millions of times 
larger than the stellar radius. Using data from PSR J0737 — 3039A/B we dramatically constrain 
the non-dynamical Chern-Simons coupling to Mcs '.= > 33 meV, approximately a hundred 

billion times better than current Solar System tests. 

PACS numbers: 04.80.Cc,04.60.Cf,04.60.Bc,04.50.Kd 



I. INTRODUCTION 

String theory is an intricate web of mathematically 
beautiful hypothesis that promises to unify all forces 
of nature. General Relativity (GR) is expected to be 
its low-energy limit with possible higher-order curvature 
corrections. To date, however, string theory remains in- 
trinsically difficult to test experimentally, because these 
curvature-corrections arc believed to be perturbativcly 
Planck suppressed. Dynamical situations with large 
spacetime curvature could lead to non-linear couplings 
and enhance such curvature corrections to a constrain- 
able realm. 

One such curvature correction is the parity -violating 
Pontryagin density, which in addition to the Einstein- 
Hilbert term defines an effective theory: Chern-Simons 
(CS) modified gravity In four dimensions, this den- 
sity is a topological term that does not contribute to 
the field equations, unless its coupling is non-constant or 
promoted to a scalar field [2|. From a string theoretical 
standpoint, the Pontryagin correction is inescapable, if 
one is to have a mathematically consistent theory that is 
anomaly-free [3, HjQ- From an experimental standpoint, 
the search for the breakage of fundamental symmetries 
can provide hints that can guide theorist toward the cor- 
rect ultraviolet completion of GR. 

CS modified gravity is also motivated from the stan- 
dard model of elementary particles and from Loop Quan- 
tum Gravity (LQG). For example, in particle physics, 
we know that an asymmetry in the number of left- and 
right-handed fermions forces the fermion number current 
to be anomalous, in analogy to the triangle anomaly Q. 
This anomaly leads to the inclusion of the Pontryagin 
density in an effective fcrmionic action Similarly, 
it has been recently found that LQG also leads to CS 
modified gravity when the Barbero-Immirzi parameter 
is promoted to a pseudo-scalar field in the presence of 



fermions [1, i, 0, [IH, El • 

The signature of CS modified gravity is the enhance- 
ment of gravitational parity asymmetry, which in partic- 
ular leads to frame-dragging modifications [H, O, [l^. 
In GR, the gravitomagnetic sector of the metric couples 
to the spin and the orbital angular momentum of grav- 
itating systems, leading to corrections in their orbital 
evolution, such as precession of the orbital plane. In CS 
modified gravity, the gradient of the coupling scalar se- 
lects a preferred direction in spacetime that corrects this 
precession. Thus, observations of gravitomagnetic pre- 
cession can be used to test the validity of the effective 
theory [11,111. 

In the Solar System, this precession correction has al- 
ready been studied for an externally prescribed (non- 
dynamical) CS coupling d. Through comparisons with 
the LAGEOS and the Gravity Probe B experiment, 
bounds have been placed on the local magnitude of the 
time derivative of this field < 10^ km or its associated 
energy scale Mcs ■=!/'& ^ 10~^^ eV near Earth-^. 

From a theoretical standpoint, the effective mass scale 
for the CS term is uncertain. While it could be as large 
as the Planck scale, it is intriguing to explore the possi- 
bility that the scale is around cosmological constant scale 
A^/** ^ 1 meV. Such intrigue arises because the cosmo- 
logical constant is an example of a quantity that accord- 
ing to string-theoretic predictions could be as much as 
120 orders of magnitude larger than the observed value, 
depending on the formulation. Since CS modified grav- 
ity is also predicted by string theory, it is interesting to 
study whether its modifications arc also observable at the 
cosmological constant scale. 



^ When converting i? into an energy scale we shall employ natural 
units. Otherwise, in this paper we use geometric units. 



2 



The weakness of Solar System constraints can be qual- 
itatively understood by focusing on the ratio of the CS 
precession correction to the GR expectation. For any 
binary system, this ratio scales as {TZext/T^md)'^ , where 
TZext and TZmd are the radius of curvature of the com- 
bined system and of either compact body respectively. 
For a binary system TZext ~ a, where a is the semi-major 
axis, and TZind ~ where R is the stellar radius. In the 
Solar System, a = i?+ -|- h, where h is the height to which 
satellites can be reliably placed in orbit, while R = i?+ is 
Earth's radius. Thus, the ratio TZext /T^ind^^ ~ h/R^ <C 
1 and the CS effect is inherently small. For a binary pul- 
sar, however, TZext /TZind ~ 0(10^), which thus enhances 
the CS effect by a factor of approximately C(10^°). 

In this paper, we study non-dynamical CS modified 
gravity in the far-field, applied to gravitomagnetic pre- 
cession. We choose to work with the non-dynamical 
theory, since this has been studied in more detail (se e 

eg. 0, a Q [S M Tl, la, i£, M, 2L 21, 2^ Mm) 

and we choose the standard CS coupling scalar -d = Test, 
where Tcs is a quantity we wish to constrain with units of 
length. We recalculate the corrections to gravitomagnetic 
precession and solve the orbital perturbation equations to 
find the CS corrected, averaged rate of change of the pe- 
riastron. Using the measurement of periastron precession 
from the double binary pulsar PSR J0737-3039A/B HI], 
we place a bound on the magnitude of the time derivative 
of the CS coupling: = Tqs ^ 6 x 10^^ km or equiva- 
lently Mcs ^ 33 meV (in natural units), much stronger 
than the previous Solar System constraint. 

The remainder of this paper deals with the details of 
this calculation and it is divided as follows: Sec. |TT] de- 
fines CS modified gravity and presents the modified field 
equations of the theory; Sec. IIIII tests the non-dynamical 
framework with the double binary pulsar; Sec. IIVI con- 
cludes and points to future research. 

We shall here employ the conventions in [23|, with 
Greek letters ranging over spacetime indices, Latin let- 
ters over spatial indices only. We work exclusively 
in four spacetime dimensions, with the metric signa- 
ture (— , -|-, -|-, -|-). Round and square brackets in index 
lists denote symmetrization and anti-symmetrization re- 
spectively, namely T(^ai3) = l{Tai3 + Tpa) and T^ap] = 
\{Tai3 — Tpa)- The Einstein summation convention is 
employed unless otherwise specified and geometric units 
(G = c = 1) are used mostly throughout the paper, ex- 
cept when relating our results to those of Ref . , where 
we use natural units (h = c = I). 



where 

•Seh = 
'S'cs = 

- 

»Smat ^ 
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-gR, 

^-d *RR, 



(2) 
(3) 



d^xV^ [g^^ (V,^) (Vpd) + 2V(^)] ,(4) 

(5) 



Equation [2] is the Einstein- Hilbert action, Eq. ([3]) is the 
CS correction, Eq. ^ is the action for the scalar field and 
the last equation represents additional matter degrees of 
freedom. In these equations, ~ IGttG, a and /? are 
dimensional coupling constants, g is the determinant of 
the metric, Vq is the covariant derivative associated with 
the metric ga(3, and R is the Ricci scalar. The CS action 
depends on the Pontryagin density *RR, namely 



It JX — Itlt — fi tl a^S . 

where the dual Riemann-tensor is 
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Idea- : 



(6) 



(7) 



and e^^'^'^ the 4-dimensional Levi-Civita tensor. 

From the action we can derive the CS modified field 
equations [J, 



G„ 



-c„ 



(5Ud 



I 



(8) 
(9) 



where G,,^ is the Einstein tensor and the C-tensor G,,„ 
is defined via [29[ 

C^'" := v„e'"'3"(''V„i?'')^ -I- v^^ *R''iP^)^ . (10) 

where □ is the D'Alembertian operator, and v^r = V 
and Vaa = VctVqI? are the covariant velocity and ac- 
celeration of the scalar field. The matter stress-energy 
tensor is T^^* and the scalar-field stress-energy is T^n) 
[see eg. Eq. (67) in [111], where 



t: 



Q/3 



VaVp - ^^gapV^v" - gaf3V{-d) 



(11) 



II. CHERN-SIMONS MODIFIED GRAVITY 
AND THE FAR FIELD SOLUTION 

Let us begin by summarizing the basic equations of 
CS modified gravity (sec eg. [1^ for a pedagogical review 
or for more details). The CS action is here defined 

by 



S = S 



EH 



5*, 



CS 



(1) 



The scalar-field potential U(i?) depends on the funda- 
mental theory from which CS modified gravity derives. 
In string theory, is a moduli field with a shift symme- 
try that forces the potential to vanish. One can generally 
employ this assumption to set the potential to zero. 

One of the main ingredients of CS modified gravity is 
the CS coupling scalar i9 = ■d(x'^), which is a function of 
spacetime. If this field were a constant, then the modified 
theory reduces to GR because the Pontryagin density is 
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a purely topological term. The coupling constant a de- 
termines the coupling strength of the CS scalar and the 
Riemann curvature, while the coupling constant /3 de- 
termines the magnitude of the energy in the scalar field. 
The choice of units for either (a, /3, ■&) fixes the units for 
the rest. For example, if [a] = where L is a unit of 
length and A is any real number, then \&\ = L^^ and 
[/?] = , where we required the action to be dimen- 

sionless when using natural units [h ~\). If instead we 
use geometric units (G = c = 1), then the action has 
units of I? ^ which if \{f\ — L"*, requires [a] = L?^^ and 
[/3] = 

CS modified gravity can be classified into two distinct 
sub-theories: a dynamical one and non-dynamical one. 
In the dynamical theory, /3 and a are arbitrary and the 
field equations are written in Eq. ([5]) and ([5]). In the non- 
dynamical theory, /3 = at the level of the action, and 
the field equations become 

Guy ^ C* = —T™^ , (12) 
= *RR. (13) 

The evolution equation for -d thus becomes a differential 
constraint, the so-called Pontryagin constraint, for the 
space of allowed solutions, while the scalar field is an 
externally prescribed quantity. 

In this paper, we shall choose to work with the non- 
dynamical theory^, since this has already been studied 
in detail and constraints (albeit weak) have already been 
placed on the strength of the correction 0, [l^, [ij, [l^ . 
In the non-dynamical framework, the functional form of 
the CS coupling scalar is not pre-determined. When 
CS modified theory was originally proposed [l|, a spe- 
cific choice was made, namely = Test, where Tcs is 
a constant with dimensions of length. A possible inter- 
pretation for -de is as an "arrow of time" since its asso- 
ciated embedding coordinate becomes = [tcs, 0,0,0]. 
From a mathematical standpoint, this choice is conve- 
nient, since it leaves the CS action time-translation and 
reparameterization invariant (see eg. [l[ for more details). 
From a physical standpoint, this choice is also conve- 
nient, since the Schwarzschild solution is automatically 
recovered for stationary and spherically symmetric back- 
grounds, and thus, most Solar System tests are automat- 
ically passed. The only tests that are not automatically 
passed are those that involve gravitomagnetism, such as 
LAGEOS and Gravity Probe B, and these experiments 
have been used to constraint tJc = ''cs ^ 10'^ km 
In Appendix 13 we present some informal arguments for 
why i!) = -de might be the only allowed functional for for 
the CS field in the Solar System, although a formal proof 
is still lacking. 

With such a choice of -d, one can solve the linearized 
field equations for the metric components. In Ap- 



^ Henceforth, wc choose a = k, following [J] 



pendix[Xl we show explicitly that the temporal-temporal 
and spatial-spatial sectors of the modified field equations 
are automatically satisfied, which implies that scalar 
gravitational perturbations are unaffected by the CS 
modification. The Oi field equations, however, are CS 
corrected, but they can be solved to linear order 
via (sec Appendix!^ for more details) 

where p is a matter density distribution, while Vi is its 
three-velocity, x is the Euclidean cross product and we 
have neglected any time dependance of the Newtonian 
gravitational potential. The vectorial solution presented 
here is similar to that found in [3, [lB| , except that here 
we consider generic density distributions. One can show 
that in the limit as p ^ m S^{x'^), Eq. (fTil) reduces iden- 
tically to Eq. (44) in [l3|j with the appropriate choice of 

''"cs ■ 

The gravitomagnetic potential presented in Eq. p4)) is 
similar to that found in One can show that if pvi is 
replaced with the stress-energy component appropriate 
to a homogeneous rotating sphere, then this potential re- 
duces to Eq. (54) in [3 with the appropriate choice of 
Tcs and to linear order. Care must be taken, however, 
when solving explicitly Eq. (|14p with a non-trivial den- 
sity distribution, as boundary terms might be required 
to ensure the junction conditions are satisfied 



III. ASTROPHYSICAL TESTS OF CS 
MODIFIED GRAVITY 

A. Weak-Field Tests 

The lack of a CS correction to the scalar sector of the 
gravitational perturbations implies that most astrophysi- 
cal process are unaffected. For example, the equations of 
structure formation remain untouched because the Pois- 
son equation is not CS corrected and the stress-energy 
tensor remains locally conserved. The vectorial sector of 
the metric, however, is CS modified in a normal direc- 
tion relative to the GR prediction. For randomly ori- 
ented velocities, the average value of the leading order 
CS correction (that shown in Eq. (|14p) in fact identically 
vanishes, simply because the correction is odd in and 
< >= 0. In many astrophysical scenarios, however, 
the velocity field is not randomly oriented. One such 
case are binary systems, where the CS correction leads 
to an anomalous frame-dragging effect. 

Anomalous frame-dragging induces modifications on a 
variety of astrophysical processes, such as the formation 
of accretion discs around protoplanetary systems and the 
evolution of neutron star spins. The CS correction, how- 
ever, would be hard to detect in such processes because 
it scales inversely with the radius of curvature of the sys- 
tem, as one can see from Eq. ((Ti)) . Galactic radii are on 
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the kpcs scale, which renders the ratio of the CS correc- 
tion to the GR prediction on the 0(10"^^) if we saturate 
the Solar System constraint 

Although the CS correction is insignificant in the evo- 
lution of non-compact astrophysical sources, this is not 
the case for binary systems. For example, inspiraling 
black hole binaries would be ideal laboratories to test 
CS modified gravity. Such systems do not radiate elec- 
tromagnetically unless surrounded by an accretion disk, 
but gravitational wave observations with space-based or 
earth-based detectors could be used to test the modified 
theory [13, HU, [13, [U ■ Such observations would be sen- 
sitive to the integrated history of the CS term, instead 
of its instantaneous value. 

Another type of binary system that shall be used in 
this paper to test the modified theory are binary pulsars. 
In such systems, there are two important scales the CS 
correction could couple to: the radius of curvature of the 
system, which is proportional to the semi-major axis a; 
and the radius of curvature of either component, which 
is proportional to the radius of either body R. As we 
shall see, the GR prediction for the precession of the pe- 
riastron scales as a"^, while the CS correction scales as 
a~^R~^, which implies that observed binary systems are 
preferred laboratories to test the modified theory since 
a/R^O{10^). 

One might worry at this juncture that the CS correc- 
tion dominates over the GR solution when M/a <C 1, 
since the former seems to decay more slowly with semi- 
major axis than the latter. In this paper, however, 
we shall assume that the solution in Eq. (fT4|) applies 
when the CS correction is small relative to the GR so- 
lution. This is indeed the case provided Tqs R^ /a?^ 
which for the binary system under consideration becomes 
Tcs ^ 10~^ meters. We shall see that the bound derived 
in this paper forces Tcs to be well below this value, and 
thus, the approximation made arc indeed valid. 



B. Binary Pulsar Test 

Consider binary systems of spinning neutron stars, 
whose orbital evolution we shall model through a 
geodesic study of a compact object in the background 
of a rotating, homogeneous sphere. Following [2|, the 
stress-energy tensor of this sphere will be described by 

To; = —ji, where the current ji = pf^ {vi x -F^ Q{R — r) , 

with po some constant density, $7' = [0, 0, f2] the angular 
velocity vector in Cartesian coordinates, = [x, y, z] the 
distance from the center of the sphere to a field point and 
R the radius of the sphere. The total mass of this sphere 
is M = 3po/(47ri?^), while its spin angular momentum 
J* = in^, where / = (2/5)Mi?^ is the moment of inertia 
The motion of the compact body in this background 
is governed by the geodesic equations a ~ — 4v x B, 
where we have neglected time-dependent scalar poten- 
tials and where a and v are the three-acceleration and 



three-velocity of the compact object. The gravitomag- 
netic field is B := W x A, where the gravitomagnetic 
potential is Ai := —gQi/A. As anticipated in the previ- 
ous section, both the field and the potential have been 
computed for this stress-energy to arbitrary order in 
Tcs [2| , but we shall here work only to leading order Tcs , 
since as we shall see second-order effects will be negli- 
gible. The gravitomagnetic field can then be written as 
Bcs = B — i?GR , with 



Be 



Co 



cos[C(r)] 



J 



tanC [J xfj - i^J- 

(15) 

whcrej(r)^= 2r/Tcs, Cq = 15Tcs/(4i?) sin[e(i?)], r = 
r/r, J = J I R^ and • the Euchdean dot product. Note 
that E g. (p5l) accounts both for the homogeneous solution 
of [3, Il5j and the boundary term found in [3| . 

From the gravitomagnetic field, we can straightfor- 
wardly compute the CS correction to the geodesic ac- 
celeration by taking the cross-product with the velocity 
vector 

Acs = -'^\cos[^{r)](^vx 

- sin[^(r)] J [v ■ f) — r {^J ■ V 

- cos[C(r)] {j ■ (iT X f)| , 



(16) 



plus subdominant terms that scale with higher powers of 
the derivatives of the CS scalar. 

One must be careful when expanding solutions in 
Tcs, since this quantity has units of length, and thus, 
corrections will arise as combinations of 0{tcs/R) and 
©(tcs/o). The approximations made so far hold pro- 
vided these combinations arc much smaller than unity. 
As we already argued, Tcs ^ R^/a^, which supersedes 
the above requirements. Care must be taken, however, 
since the argument of the oscillatory functions scales as 
1/tcs, and thus, any spatial derivatives of Ai (or goi) will 
be larger than Ai (or goi) by one power of Tcs- When 
this fact is taken into account, the results for the grav- 
itomagnetic field found in [3, [111 ^'^'^ are in formal 
agreement'^. 

We shall here parameterize the trajectory of the com- 
pact object in terms of equatorial coordinates. We shall 
thus define the triad: 



r — [cosu, costsmu, smism uj , 
t ~ [— sinu, cos i cos u, sin t cos u] 
ri = [0, — sint, cosi] , 



(17) 



to describe radial, transverse and normal directions rel- 
ative to the comoving frame in the orbital plane. The 



^ As already discussed, 
is not modeled in Il4l 



131 obtains a boundary contribution that 
1511 . because the latter employed a point- 
particle approximation, while the former dealt with extended 
bodies. 
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quantity t is the inclination angle, w is the argument 
of periastron, u = f + w with / the true anomaly and 
51 = is the right ascension of the ascending node, cho- 
sen in this way so that the line of nodes is co-aligned with 
the X vector |32| . 

The perturbation equations for the variation of the Ke- 
plcrian orbital elements is governed by the projection of 
the geodesic acceleration on this triad. To leading order 
in Tcs, however, only the radial projection := a • f is 
CS modified, leading to a^^ := — a^'^: 

a^^ = —AcquJ' {cos L cos (r)] + sin l cos u sin (r)] } . 

(18) 

We can now compute the variation of the orbital ele- 
ments by studying the perturbation equations (s^ 



da 
'dt 
de 
H 
di 
di 

dn 

m 

d(jj 
Itt 



nVT 



e ar sin f + at — 
r. 



ar sin f + at 



1 



cos / + - 1 



avT" 



=a„— cosM, 
°2 a 



1 



=a„- smu, 
=^ a 



n a sin — 

-ar cos / -t- at ( 1 



= n —ttr'- \/l — f ^ + cos ] . (19) 



n a e 

dVL 
OS i — 
dt 

2 



n a 



duj 



- sin / 



dVl 



dt 



dt 



where n = Mja? is the unperturbed Keplerian mean 
motion, e is the eccentricity, and M. is the mean anomaly. 



Since a<=" • t = ©(r^J 



n. fir 



while Wcs 



—arjin a e) cos /, to leading order in the eccentricity. 

The average of the rate of change of w can be computed 
by integrating cj over one orbital period: 



< w >: 



-dt 



27r 



1 



27r(l + ecos/) 



r#, (20) 



during which we shall assume the pericenter 
is approximately constant, so that ii ^ f = 



-3/2 



and the motion of the 



7i(l + ecos/)^ (1 - e^) 
compact object can be described by a Keplerian el- 
lipse, where r = a (l — e^) (1 -I- ecos/)~^. This last 
assumption is justified by the fact that the motion of 
test particles about any arbitrary background remains 
unchanged relative to the GR prediction, ie. test-particle 
motion satisfies the geodesic equation both in GR and 
in CS modified gravity. Such is the case provided the 
strong-equivalence principle is satisfied, which is guaran- 
teed by the Pontryagin constraint in the non-dynamical 
version of the theory or by the scalar field equation of 
motion in the dynamical version (see [sij for a proof). 
Finally, the integrals in Eq. ([20|) shall be approximated 
with a small eccentricity expansion e ^ 1. 



The averaged rate of change of the periastron can then 
be decomposed into a GR prediction plus a CS correction, 
where the latter is given by 

15 J Tcs ^ . f2R\ . (2a\ 

A sm sm , (21) 

V'^cs/ \TcsJ 



and where the projected semi-major axis is X := asint. 
Equation (pi]) neglects terms of order unity, since the 
dominant contribution scales as and we here concen- 
trate on systems with small but non- vanishing eccentric- 
ity. The scaling occurs because although w scales 
as cos/, so does aj;^, and thus, the leading order term 
in e does not vanish upon integration, unlike in the GR 
case. The orbital orientation, however, is ill-defined for 
exactly zero eccentricity, and thus, the limit e ^ is 
meaningless. The scaling in the precession of the pe- 
riastron of Eq. (PT|) is consistent with other precession 
results studied in the Solar System As discussed in 
the Introduction, the ratio of CS correction to the GR 
expectation scales as a'^Tcs/R^, since (w)^,^^ ~ J/a^. In 
the Solar System, however, a/R is very close to unity, 
while for binaries a/R - O(IO^). 

At this junction, one might worry that the calcula- 
tion of the CS modification to w is not sufficient to place 
a bound on the non-dynamical theory with the binary 
pulsar, since other relevant quantities that play an im- 
portant role in the test might also be CS modified [l^. 
Although this is generally true, the CS modification to 
other quantities turns out to be either identically zero 
or subleading. This is so because Cqo identically van- 
ish for 7? cx t in the non-dynamical formalism, while 
Cij is proportional to gravitational wave perturbations 
only. These waves will be CS modified, but the CS cor- 
rection to the quadrupole formula is subleading [hf-^ cx 
{tqs/ DL)d^Iij / dt^ , where hj is the quadrupole moment 
and -Dl is the luminosity distance to the source (sij]. 
and thus, the CS correction to the rate of change of the 
binary period P is also subleading. Since the only post- 
Keplerian parameters that is CS corrected to leading or- 
der is tj, we can treat, as a first approximation, all other 
orbital quantities as given in 

With these remarks in mind, observations of the pre- 
cession of the periastron in the double binary pulsar 
PSR J0737 - 3039A/B [H can be used to test CS 
modified gravity. We shall treat pulsar A as the ro- 
tating homogeneous sphere, and pulsar B as the test 
body in orbit around the sphere, where the bodies are 
sufficiently separated that we can neglect tidal interac- 
tions. The relevant system parameters are [s^ the mass 
Ma = M k, 1.34M0, the projected semi-major axis 
X w 1.41 s, the eccentricity k, 0.088 and the incli- 
nation angle l « 89(— 76, -1-50) deg. From the projected 
semi-major axis we can deduce that af, 4.24 x 10^ km, 
where we used the nominal value for the inclination an- 
gle. We assume the standard value for the moment of 
inertia of body A / w 10'^® kg m^, which leads to a ra- 
dius of Ra w 9.69 km and an angular momentum of 
J A ~ 2.8 X 10^° kg m^ s~^, since the rotational pc- 
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de/dt [km] 

FIG. 1: CS correction to the precession of the periastron as 
a function of res for the system parameters of PSR J0737 — 
3039A/B. 

riod has been determined to be 22 milliseconds [s^. IsH- 
The precession of the periastron has been measured to 
be ti; ~ 16.96 degrees per year, in complete agreement 
with the GR expectation, with an overall uncertainty of 
approximately S = 0.05 degrees per year"' [s^. [35|. 

We can then constraint Tcs by requiring graphically 
that {u!cs) be less than 6. The uncertainty in geomet- 
ric units becomes S = 9.2 x 10~^'^radians km~^. Fig- 
ure [1] presents a plot of (wcs) as a function of Tqs- A la 
constraint can then be derived from this figure, namely 
^ = Tcs < 6 X 10"^ km, or simply Mcs := = t'^^ > 
33 meV in natural units, which is approximately 10^^ 
times stronger than current Solar System constraints. 

We have checked that terms higher order in e or Tcs do 
not significantly affect this bound, which is however pri- 
marily affected by uncertainties in the semi-major axis. 
Even with the most pessimistic choice of a, the bound 
deteriorates only by a factor of twenty, still leading to a 
constraint 10^° times stronger than the Solar System one. 
Also note that this bound is consistent with the approx- 
imation made to derive the CS correction to periastron 
precession. 



IV. CONCLUSIONS 

We have studied non-dynamical CS modified gravity 
with a time-like CS coupling scalar. Until now, the 
only constraint on non-dynamical CS modified gravity 
{Mcs ^ 10~^^ eV) came from Solar System experiments 
due to CS corrections to frame-dragging 0, [ij, [l^ . We 
here calculated the leading-order CS correction to post- 
Keplerian parameters of binary systems. We find that 



* The precession of the periastron has been measured to higher 
accuracy for pulsar A, but we adopt here the larger uncertainty 
so as to derive a conservative bound on mcs- 



the precession of the periastron is the only parameter 
that is CS corrected to leading order. This corrections is 
such that its ratio to the GR expectation scales as a^/i?^, 
where a is the semi-major axis and R is the neutron star 
radius. For the binary pulsar considered here, this ratio 
is of 0(10^") which leads to an enormous enhancement 
over previous Solar System constraints: Mcs > 33 meV. 
This constraint is approximately a hundred billion times 
stronger than current Solar System constraints. 

Although this paper constrains the non-dynamical 
framework of CS modified gravity to unprecedented lev- 
els, it cannot do the same for the dynamical formulation. 
Corrections to post-Keplerian parameters in the dynami- 
cal theory are of high post-Newtonian order, because the 
Pontryagin density vanishes to leading order. Meaning- 
ful tests of the dynamical formulation would then have 
to rely on strongly gravitating sources. 

One such scenario is the inspiral and merger of compact 
objects. Dynamical CS modified gravity should correct 
both the trajectories of such objects as well as the gener- 
ation of gravitational waves. A detection of such waves 
with LIGO or LISA could then be used to plac e st ring ent 
bounds on the dynamical formulation (20l. l2ll. [sol . l3l| . A 
program that pursues just such a calculation is currently 
ongoing [10, [3l| . 
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APPENDIX A: FAR-FIELD SOLUTIONS IN 
NON-DYNAMICAL CS GRAVITY 

In this Appendix, we study solutions to the modified 
field equations in the far-field. We begin by describing 
the approximation schemes employed and then tackle the 
modified field equations order by order. We shall here ini- 
tially allow to be arbitrary, but we shall present some 
informal arguments that suggest might be the only so- 
lution allowed in the Solar System, albeit a formal proof 
is lacking. 

1. Approximation Schemes 

Consider the far-field expansion of the line clement 

ds^ = -{l + 2(j))dt'^ + 2gidtdx' + {l-2'i{;)S,jdx'dx^ , (Al) 

where t, a;* are Cartesian coordinates, 6ij is the Euclidean 
metric, {ip, (j)) and gi := goi are scalar and vectorial 
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perturbation potentials respectively in the longitudinal 
gauge = 0). 

The perturbation potentials and the matter sources 
that generate them shall be treated perturbatively, in 
a post-Newtonian (PN) sense, where the latter are as- 
sumed slowly- moving (e := v j c <^ 1), weakly-gravitating 
and isolated (sec eg. [Sj])- We shall ignore gravitational 
wave perturbations, since these have been partially stud- 
ied elsewhere [ll, H S ISi. In particular, we shall 
model these sources via a perfect fluid stress energy ten- 
sor, such that Too = P = 0{e'^), Toi = -pVi = 0{v^) 
and Tij = 0{e^), where p is density and is the three- 
velocity of the fluid. 

We shall be concerned here with binary systems, whose 
exterior gravitational field (or metric) can be modeled in 
GR (to leading post-Newtonian order) via 

mi m2 
i>GR ^ 0GR ^ 1 , 

g'o. - —vm + '^vm, (A2) 

where mi_2 are the component masses, ri^2 = \x^—x\ 2(i)| 
are their field point positions, with x\ 2{t) the trajecto- 
ries and v\ 2(<) = x\ 2{t) the velocities. Note that these 
potentials are not spherically- symmetric and are here ex- 
pressed in Cartesian coordinates, in spite of the appear- 
ance of fiducial radial distances ri_2- Such potentials be- 
come spherically symmetric only in the limit TO2 — > 
and mi fixed (or viceversa), in which case one recovers a 
boosted, Schwarzschild metric in harmonic coordinates. 
If these objects are spinning, then the gravitomagnetic 
sector of the metric g^.^ acquires more terms prop ortional 
to the spin angular momentum (see eg. p^. Il5l|). These 
equations can be derived by assuming a point-particle 
approximation, but we relegate any such details to the 
post-Newtonian reviews in [sl, [sj . 

We shall concentrate on a rather special, yet physi- 
cally reasonable subset of metric perturbations: poten- 
tials that are small CS deformations of GR solutions. 
In other words, on top of the PN perturbative scheme, 
we shall employ a small- coupling approximation. In the 
dynamical formalism, this can be achieved by expand- 
ing in C := /{k(3M^) <^ 1, where M is the gravita- 
tional mass (a length-scale) associated with p. In the 
non-dynamical scheme, the perturbation parameter be- 
comes C = \d^-d\/M <C 1, which for i? = i?c becomes 
§/M ^ Tcs/M <C 1. 

The combined use of a PN expansion and the small- 
coupling approximation defines a bivariatc perturbation 
scheme, where both e and C, can be treated as indepen- 
dently small parameters. Moreover, in the dynamical 
framework, this scheme defines a boot-strapping frame- 
work in which one can first solve the evolution equation 
for -d in the non-CS corrected background, and then use 
this •& to solve the modified field equations to first-order 
in the CS correction. For more details on this boot- 
strapping scheme or the small-coupling approximation 



as applied to the dynamical theory, we refer the reader 
to [sq. 

Based on these considerations, we shall commonly 
make the decomposition A = A'^^-\-X(a)^ where A is any 
metric perturbation, A'-^^ is a GR solution of ©(C") and 
XA is some undetermined potential of 0{C,). Moreover, 
we shall require that A and XA be at least of the same 
order in e, such that we can search for CS-deformed solu- 
tions. The metric perturbations shall then be expanded 
as 

^ = V'gR + X(4>) , (t> = <l>GR+ XW > = 5gR + 

(A3) 

where ^pcR and ^g/? are both of 0{e'^X^); 5gr of 
C(e^,C°), and X{4>) are both of O(e^0, and x|g) 

is of 0{e^, C). The notation 0{e"^,C) stands for a term 
of 0(e™) or C(C"). 

Such a decomposition neglects CS corrections that 
modify the leading-order behavior of GR. For exam- 
ple, we shall not consider a perturbation X{4>) of 0{e, ^). 
Such an assumption is justified on the basis of the small- 
coupling approximation and the fact that GR has been 
found to agree to incredible precision with Solar System 
experiments (33| . The above expansion then guarantees 
that to lowest order all such experiments are passed, with 
the new potentials X()/'),(0),(s) leading only to small per- 
turbations. 



2. Scalar Metric Perturbations 

The scalar metric perturbations are traditionally 
solved for by studying the field equations of a theory to 
0{e'^). For an arbitrary -d, the modified field equations 
become 

= 47rp, (A4) 
l^'' [(d,^) dkV^ + (9,„i?) d^k] (V' + = 0, (A5) 

(%V2 - d,j) (0 - V^) = - (did) e(>'9,)fc {i: + 0)(A6) 

where is the Laplacian in Cartesian coordinates, 
gufc ._ gOij/c g^-^i^ pjk jg ^j^g Levi-Civita symbol. 

The decomposition of Eq. (jASp simplifies the modified 
field equations. To 0(e^,^"), the field equations reduce 
to those of GR, namely 

VVgr = 4^P, (A7) 
(5y - (</,oK - Vgr) = 0, (A8) 

where the temporal-spatial component automatically 
vanishes. E quat ion (|A7[) leads to the standard GR so- 
lution [3,[il,lil 

V'GROc-|^^dV, (A9) 
while Eq. (jASP implies 0gr = V'gr- 
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To next order [0{e^X)i 
become 



V'X(^) = 0, 



the modified field equations 



~ jk 

^7 



0, 



Once more, when we apply the CS deformed decom- 
position of Eq. (|A3p . the modified field equations sim- 
plify. To 0{e^, C"), only the temporal-spatial component 
(AlO^f these equations survive: 

(All) 



(A17) 



~ kl 



d.)ki^c 



Equation (|A10[) forces X{if!) to be a solution to the Laplace 
equation. The remaining two equations [Eqs. (|A11[) 
and (|A12[) ] are clearly satisfied if = dc- 

We have searched for other choices for that satisfy 
these equations and we have found the following sufficient 
conditions: 



(Al^hich leads to the exterior gravitomagnetic solution of 
Eq. ^ (see eg. BEIH). 

To next order [©(e'^,^)], the field equations become 



(i) a,^? = 0, (ii) d,^ = 0{e). 



(A13) 



Clearly, the choice d = -de satisfies either (i) or (ii). Op- 
tion (i) forces both terms in Eq. (jAlip to exactly vanish, 
as well as the right-hand side of Eq. (jA12[) . while op- 
tion (ii) forces these terms to vanish only perturbatively, 
since then (9i?)(9'i/'GR) = 0{e^,C)- Either option then 
forces X(4:) = X(0): but since these functions must satisfy 
the Laplace equation and the metric (and thus X(i>},{4>}) 
must be asymptotically flat at spatial infinity^ , wc choose 

XW =O = X(0)- 

Although the conditions presented above are sufficient, 
we cannot formally prove that they are necessary to sat- 
isfy the modified field equations. In other words, al- 
though we have failed to find a solution to the above 
system of differential equations, we cannot prove that a 
solution does not necessarily exist. 



3. Vectorial Metric Perturbations 

The vectorial sector of the metric perturbations can be 
solved for by studying the field equations to 0{e^). For 
an arbitrary scalar field, these equations become 



(9,1?) v^a.g,. = (du^) di.ok, 

4 (did) 



(A14) 



~ijk 



did 



{did) h 



{d,,gk) , (A15) 

(a„fcz9)9,[,)gfc], (A16) 
where we have used the longitudinal gauge condition. 



- kl 



n~ nl , 



Since the metric is asymptotically flat, it must decay as 1/ri 
or l/r2 at spatial infinity, but the Laplacian of these functions 
is non- vanishing. For example, in the limit of vanishing 7712, 
the Laplacian of 1/r becomes proportional to the Dirac delta 
function. 



- kl 



dh^'V^dugr 
~ijk 



(A18) 



[did) {d.^gD , (A19) 

(did) e~(/^ {\/^dkg,)ar.) 

26(,"'(9„,^)5,[,)gf , (A20) 

With the choice d = dc, not all the field equations are 
automatically satisfied, with the temporal-spatial com- 
ponent becoming 



(9,^) {d,,gn + ~^' 
4 (did) e(,'^- (a.-jfeV-GR 

^ (dnkgr 



~ kl 



^ kl 9 

2 



dkg'^ 



(A21) 



whose solution is given in Eq. (|14p by noting that t/oi = 
gi. If a non-trivial density distribution is used, such as 
the homogeneous sphere in Sec. IIIIl then one must en- 
sure that the solution to Eq. (|A21[) accounts for possible 
boundary contributions that arise to guarantee the junc- 
tion conditions are satisfied. 

Let us now argue that the conditions of Eq. (jA13P 
together with the modified field equations to this order 
lead directly \,o d =- d c- With either condition, one can 
show that Eq. (|A18|) and the second hue of Eq. (|A19|) 
either automatically vanish or become of C(e^, Q. Equa- 
tion (|A20p . on the other hand, leads to li^i'^'d {dj)kg?^) = 
0, which forces d = since gf^ does not vanish. Com- 
bining Eq. (KT3\i with d ^0 one is then led to d = dc- 

Once more, although the conditions in Eq. (|A13p are 
sufficient, we have not succeeded in formally proving that 
they are necessary to satisfy the modified field equations. 
In other words, we could not mathematically prove that 
Eqs. (|A18[) - (|A20[) does not possess some other obscure 
solution that we have missed. 



4. The Pontryagin Constraint 

One last issue to consider is whether the solution found 
above satisfies the Pontryagin constraint of the non- 
dynamical theory. The Pontryagin density is indepen- 
dent of the CS scalar, but for the line element in Eq. (|Aip . 
it is given by 

*RR^ -4?'-'"'= (d/g,) dki (</> + V') + O(e'). (A22) 



With the GR solutions of Eq. (|A2p , the Pontryagin den- 
sity identically vanishes to this order, since diri^2 = '^12 
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and the Levi-Civita symbol is completely antisymmetric 
under index permutation. If one were to include spin 
correction to the gravitomagnetic components, then the 
Pontryagin density would not vanish. In this sense, the 



Pontryagin density, to leading order, is of O(e^), and 
thus, unless a post-Newtonian expansion is carried to 
high order, this constraint docs not affect the analysis 
of this paper. 
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